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All-dielectric metamaterials consisting of high-dielectric inclusions in a low-dielectric matrix are 
considered as a low-loss alternative to resonant metal-based metamaterials. In this contribution we 
investigate the applicability of the point electric and magnetic dipoles approximation to dielectric 
meta-atoms on the example of a dielectric ring metamaterial. Despite the large electrical size of high- 
dielectric meta-atoms, the dipole approximation allows for accurate prediction of the metamaterials 
properties for the rings with diameters up to 0.8 of the lattice constant. The results provide 
important guidelines for design and optimization of all-dielectric metamaterials. 


I. INTRODUCTION 


In recent years metamaterials have aroused great at¬ 
tention as a tool for achieving desired electromagnetic 
properties, which might not exist in nature Qi- Among 
them a special attention has been paid to all-dielectric 
metamaterials [bHil , which do not suffer from the Ohmic 
losses contrary to the metal-based ones. All-dielectric 
metamaterials with various constitutive units, such as 
spheres [9|, cubes 0 and rods [nllil , were shown to 
exhibit the negative magentic permeability as well as neg¬ 
ative pIMTbj l and zero refractive index [iQ. Silicon nan¬ 
odimers, exhibiting both electric and magnetic dipole re¬ 
sponses, have recently been experimentally demonstrated 
to enhance electric and magnetic fields in the optical 
range [13. High permittivity dielectric rings are of par¬ 
ticular interest, because they possess a strong broadband 
magnetic response [iiil. In addition, they have more 
degrees of freedom for tuning the resonance frequency by 
geometrical parameters than, for example, spheres, cubes 
or cylinders. Furthermore, they exhibit resonances on 
the THz frequencies |19lj^ . a very promising range for 
emerging applications. 

Being as the most frequent appearence periodical 
structures composed of subwavelength elements (the 
metamaterial lattice constant is on the order of A/20 — 
A/4), three dimensional metamaterials are often treated 
as continuous media (the procedure referred to as ho¬ 
mogenization), described by some effective parameters, 
e.g. permittivity s, permeability /r, refractive index n 
and impedance 2 :, that simplifies their description. Var¬ 
ious methods of homogenization have been suggested 
based on, for instance, field averaging, single-interface 
scattering, non-local dielectric function and other ap¬ 
proaches (see the brief overview of the methods in @). 
The most popular practical way of retrieving the effec¬ 
tive parameters remains the Nicolson-Ross-Weir (NRW) 
or S-parameters method based on inversion of formulas 
for Fresnel’s reflection and transmission coefficients. It 
was originally proposed for conventional media and 
then extended to metamaterials [25| . The NRW method 


was accused of the loss of the retrieved effective param¬ 
eters’ physical meaning, since they do not satisfy the 
locality principle 
’’anti-resonance” 


^6H2c 


and presence of the so-called 
(non-physical for passive materials 
lm{e) < 0 and Im(/i) < 0; we assume optical exponen¬ 
tial convention exp(—icjt)). Moreover, introduction of 
bulk effective parameters to metasurfaces or few layers 
metamaterials is questionable. Obviously the properties 
of a single layer or a layer with only few neighbors can 
be substantially different from the properties of a layer 
with many neighbors. In some cases of strong interaction 
between meta-atoms it is not possible to introduce any 
meaningful effective parameters [3l|. 

We use another approach in which meta-atoms are re¬ 
placed with point electric and magnetic dipoles 113 Hi, 
[ 33 ]^ . The advantage of such approach is that the in¬ 
teraction of dipoles can be explicitly taken into account 
either numerically for the arrays up to few hundreds by 
few hundreds elements or analytically for an infinite array 


However, the applicability of the point dipoles approx¬ 
imation to all-dielectric metamaterials has not been ver¬ 
ified yet. In order to achieve strong resonant features, 
high refractive index dielectrics are typically used, and 
thus even geometrically subwavelength meta-atoms may 
become electrically (optically) larger than a wavelength 
exhibiting rapid field variations within their volume. 

The goal of this work is to verify the applicability of 
the electric and magnetic dipole approximation for de¬ 
scription of all-dielectric metamaterials. We make the 
study on the example of dielectric ring metamaterials. 
We start with considering a single layer of meta-atoms 
and investigate the dependence of their electric and mag¬ 
netic polarizabilities of the geometrical parameters. Then 
in order to figure out the limits of the dipole approxima¬ 
tion we show how the retrieved polarizabilities depend on 
the metamaterial lattice constant. Finally we consider a 
few layers of dielectric rings and compare the analytical 
predictions for reflectance R with the full-wave numerical 
simulations. 

The paper is organized as follows. Sec. [II] is devoted 
to research methodology. The results of numerical sim- 
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ulations and its comparison with analytical predictions 
are presented in Sec. IIIIl Discussion and conclusions are 
assembled in Sec. m 


II. METHODS 


We consider dielectric rings of the outer and inner radii 
R 2 and Ri , respectively, and thickness t arranged in the 
square lattice with lattice constant d within each layer 
and distance between the layers h [see the geometrical 
parameters in Fig. [Ha)]. The rings consist of dielec¬ 
tric with material permittivity 5 = 100 + i suspended in 
free space. A plane wave is normally incident along the 
2;—axis and its electric field directed along the x—axis. 
We are interested in the lowest order resonances of the 
rings classified as magnetic and electric dipole modes [see 
Fig. [T][b) and (c)]. 

In order to obtain electric and magnetic dipole polar¬ 
izabilities de and am for an individual meta-atom, we 
treat every ring in the array as a combination of point 
electric and magnetic dipoles [32| oriented along the x— 
and y— axes, respectively, in correspondence with the 
polarization of the incident wave. The equivalent dipoles 
are positioned in the center of the rings. The electric 
p and magnetic m dipole moments are induced by the 
external field (amplitude of the electric and magnetic 
field are Eq an d i^o, respectively and Eq = rjoHo, where 
Vo = v/W^ ~ 1207r[f]] is the free-space impedance) 
as well as by all the neighbors within the infinite array 
(interaction field Einter ^ Hinter)- In the case of a single 
layer, taking into account the identity of the dipoles, we 
express p and rh 

P — a^i^Eo “h Einter) — de(-£/0 “1“ ^od [3 p), (1) 

rh = am{Ho + Hinter) = am{Ho + iJ.o(f/3°rh), (2) 

where the interaction constant is 
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and i?o=d/1.438 [33l. 

The amplitude reflection and transmission coefficients 
can be split into electric Re and magnetic Rm terms js^: 
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The polarizabilities of the particles are then expressed as 
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We introduce the normalized dimensionless polariz¬ 
abilities ae = de/^od^, am = am/Vod^. Then the inverse 
polarizabilities are 


J_ = «0 , ikd J_ ^ oO _ ikd 
ae ^ R + T-1' am ^ R-T + l 


( 7 ) 


Considering two layers of particles one should take 
into account the interaction between the layers. The 
wave propagates as shown in Fig. 1(a), and the coordi¬ 
nate system origin is placed in the center of a ring in 
the first layer. The propagation occurs towards nega¬ 
tive z, E{z) = Eoexp{—ikz). The system of equations 
for the normalized dimensionless dipole moments 
Pi ,2 = Pi, 2 /{£od^Eo),rni ^2 = ^ 11 ^ 2 /{pod^Ho) reads 

—pi = 1 + /3°pi + /3“'*P2 + I3~mrn2, (8) 
ae 

—P 2 = exp{ikh) + 13° P 2 + P^pi + ( 9 ) 

ae 

-TOl = 1 +/3 °TOi +/3“^m2 +/3“eP2, (10) 

—m 2 = expiikh) + 13°m 2 + /3'*mi + PmePi, (H) 

^m 

where the interaction coefficients can be found with ap- 
proximate formulas (valid up to kd < 1.5) 
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and the following equalities are valid 
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Thus for two meta-layers we get the system of linear 
equations for unknown dipole moments: 
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After finding the dipole moments one can calculate 
the electric and magnetic contributions to reflection and 
transmission coefficients 

Re = ^ikd{pi+p2e"^R, (15) 

Rm = -^ikd{mi + m2e®*^'*), (16) 

Te = ^ikdipie'-'^^ + P 2 ), (17) 

Tm = + m2), (18) 

and finally obtain R and T 

R = ^ikd{pi - mi - m2e*^^), (19) 

T = + ]^ikd{pie^^^ +P 2 + + m 2 ). (20) 































3 


For larger number N of layers the procedure is analog¬ 


ical. At first one finds the dipole moments solving the 
system of equations 


ph ^0 _ X 

l^2h ph 

_niN-2)h _niN-3)h 
Hem Hern 
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Hem Hem 
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and then calculates the reflection and transmission coef¬ 
ficients 


i?= (22) 

n=l 

1 ^ 

T = + -ikd^^iPn + (23) 

We performed numerical full-wave simulations with 
CST Microwave Studio [35| in the time domain with ef¬ 
fectively periodic (perfect electric along x— and perfect 
magnetic along axis) boundary conditions and nor¬ 
mally incident plane wave excitation. 


III. RESULTS 

For the rest of the paper we consider the thickness of 
the rings fixed to t/d = 0.2, distance between the centers 
of the rings in neighbor layers h = d and variable inner 
and outer radii Ri, ^2- The frequency range is 0-0.3 d/X 
in order to stay within the metamaterial regime. To be 
noted, formulas m,m for the interaction coefficients 
are applicable for kd < 1.5 which corresponds to 
d/X < 0.25. 


o{N-l)h 

Hem 


oiN-2)h 

Hem 

oiN-3)h 

Hem 


Pi 

P2 

P3 


I 

^ikh 

^i2kh 

ph 


rriN-i 

_ rriN _ 


^ikh(N-2) 

^ikh(N-l) 


are shown in Fig. mb) and (c), respectively. Conse¬ 
quently, the retrieved electric [Fig. [2][b)] and magnetic 
[Fig. [2][c)] polarizabilities exhibit resonance behavior at 
the frequencies of the reflection maxima. Increase of the 
size of the rings expectedly leads to decrease in the res¬ 
onant frequency and increase of the resonance strength 
[Fig|2l[b),(c)]. Apart from the low-order dipolar modes 
some higher order modes can be excited (see shaded area 
in Fig. [2]), however, they are lying in the high frequency 
zone and will be excluded from the further analysis. 

When R 2 is fixed, but inner radius Ri is varied, one of 
the reflection peaks stays almost unchanged whereas an¬ 
other one shifts dramatically [Fig. [3][a)]. Increasing Ri 
influences strongly on the location of the electric dipole 
peak [Fig. [3][b)], but has weaker impact on the magnetic 
dipole [Fig. [3](c)]. The magnetic field of the excited res¬ 
onance is induced by the circulation of the displacement 
current. Accordingly to the mode structure [Fig. 1(b)], it 
is obvious that the dominating contribution is given by 
the periphery of the ring. Overall the modes blue-shift 
with the inner hole radius Ri growth [Fig. [3l[b) and (c)] 
as the result of the decreasing high-dielectric material 
filling fraction. 


B. Applicability of the dipole approximation and 
variation of the array period 


A. Variation of the outer R 2 and inner Ri radii 

In this section we investigate the influence of geomet¬ 
rical parameters, namely, the outer and inner radii on 
the resonant properties of the dielectric ring metamate¬ 
rial. Let us fix the inner radius Ri/d = O.I and vary 
the outer radius in the range R^jd = 0.2 — 0.5 (Figl2]). 
The resonant peaks and dips are clearly seen in the re¬ 
flectance spectra, for example, follow the blue line in Fig. 
Ha). They correspond to the lower-order modes of the 
dielectric ring resonators, namely, the magnetic dipolar 
(lower frequency) and electric dipolar (higher frequency) 
modes. The corresponding electric field distributions 


For retrieval of ring electric olq and magnetic polar¬ 
izabilities, which are the characteristics of a single ring, 
we use formula (7) and simulate reflection and transmis¬ 
sion for an infinite periodic array of the rings. The in¬ 
fluence of the infinite number of neighbors is taken into 
account through interaction constant for which ac¬ 
curate analytical formulas exist for the square array of 
point dipoles [s^. At the same time, the rings we con¬ 
sider here are not point-like, but rather comparable to 
the period of the array. This may lead to both devia¬ 
tion of the field distribution of a certain dipole from the 
field of the dipole located in its center (thus, influence on 
polarizabilities and am) as well as on the interaction 
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field of all other dipoles (thus, influence on interaction 
constant /3^). We assume that interaction constant 
is the same as given by formula (7), and thus all effects 
related to the finite size of the rings are included into re¬ 
trieved polarizabilities ae and am- Then, retrieving the 
polarizabilities for the arrays of various periods we are 
able to characterize how different from the point dipoles 
the meta-atoms are (ideally the retrieved polarizabilities 
should not depend on the array period). We call the ac¬ 
tual period of the structure d* in order to differ it from 
the previously used constant d, while the latter is now 
used for normalization of frequency and polarizabilities 
only. 

Period d* influences the restored electric and magnetic 
polarizabilities (Fig. |4]) and they converge to a certain 
value with increasing period d*. Such behavior is natural, 
since the larger the period the smaller is the radius-to- 
period ratio i? 2 /d* and, thus, the accuracy of the dipole 
approximation improves. Generally convergence of mag¬ 
netic polarizabilities am [Fig. Ill[a’)-(d’)] is faster than 
electric ae [Fig. Ill[a)-(d)] what can mean that magnetic 
interaction of the rings is weaker than electric. In case 
of touching rings (d* = 2 R 2 ) the retrieved polarizabil¬ 
ities differ much from the asymptotic value, clearly ex¬ 
hibiting failure of the point dipoles approximation. In¬ 
deed, touching rings can hardly be considered as a set 
of separate rings, but rather as a periodic subwavelength 
grating or a perforated dielectric membrane. The polar¬ 
izabilities convergence is good for the rings radii up to 
i? 2 /d = 0.4 [Fig. |4l[a)-(c) and (a’)-(c’)]. For larger rings 
with i? 2 /d = 0.5 no convergence is observed for the pe¬ 
riod d* up to d*/d = 1.4 that lies on the boundary of 
the metamaterials regime and corresponds to d* 0.5A 
at the normalized frequency d/A = 0.3 [Fig. ^d) and 
(d’)]. Overall we may conclude that dipole approxima¬ 
tion is consistent for the rings radius up to i? 2 /d = 0.4 
and lattice constant d*/d > 0 . 8 . 


C. Several layers structure 

Knowing polarizabilities ae and am of individual rings 
retrieved from the simulated reflection and transmission 
for a single layer, we may analytically predict the re¬ 
flection and transmission for several layers of rings (see 
the methodology section for details). Direct comparison 
of these analytical predictions to numerical simulations 
can give us additional evidences of applicability or non- 
applicability of the point dipoles approximation. We use 
the same rings geometrical parameters as in the previ¬ 
ous subsection. Fig. [5] shows the reflection for different 
R 2 for 2 , 3 and 4 layers structure with Ri/d = 0.1 and 
i? 2 /d = 0.2,0.3,0.4 and 0.5. 

For smaller radii R 2 /d = 0.2 and 0.3 and even for 
i? 2 /d = 0.4 the analytical prediction is in a good corre¬ 
spondence with the numerical simulations providing an¬ 
other confirmation that the dipole approximation works 
well [Fig. [5][a)-(c)]. The discrepancy between the pre¬ 


dicted and simulated results increases with R 2 and for the 
case of touching rings R 2 /d = 0.5 [Fig. [Sfd)] the reflec¬ 
tions given with two methods do not resemble each other 
reasonably, especially in the region of resonant peaks and 
dips. Interestingly, the reflection is different not only in 
the region of polarizabilities resonances (see Fig. [ 2 ]), but 
also in other regions of reflection dips which we interpret 
as appeared due to Fabry-Perot resonances (increasing 
R 2 leads to the high-dielectric materials Ailing fraction 
increase and thus to the increase of the overall optical 
thickness of the metamaterial slab). 

We may conclude that the electric and magnetic po¬ 
larizabilities restored from a plane square-lattice array 
are able to predict the several layers reflectance for 
i? 2 /d < 0.4, but fail for larger radius. 


IV. DISCUSSION AND CONCLUSIONS 

We have analyzed electric and magnetic polarizabilities 
of the high-dielectric rings approximating them as point 
dipoles. The geometical parameters and the free-space 
wavelength are normalized to the unit cell size, thus the 
results of our analysis are wavelength-independent and 
virtually are applicable to any freqeuncy range. 

Variations of geometrical parameters allow for tuning 
the resonance frequencies for ae and am nearly indepen¬ 
dently: outer radius R 2 influences the resonance frequen¬ 
cies of both electric and magnetic polarizabilities (Fig. [2]) 
while the inner radius predominantly influences electric 
polarizability ae (Fig. [3]). 

We come to a natural conclusion that decreasing the 
structural unit size leads to a better prediction of the elec¬ 
tromagnetic properties with the point dipoles approxima¬ 
tion. The dipole approximation works well for the rings 
up to R 2 /d = 0.4 (Fig. 0]). This value is surprisingly 
large, since the electrical size of the ring at the frequency 
d/A = 0.15 is 2 R 2 n ~ 1.2A, thus the ring is not electri¬ 
cally small anymore and the held varies rapidly in space 
within the ring. 

For larger rings the situation is different. The rings 
with the radius R 2 /d = 0.5 are not isolated - they are 
literally touching each other. In this case the dipole ap¬ 
proximation is not valid anymore. The retrieved polariz¬ 
ability changes dramatically with the lattice constant of 
the array [Fig. |4](d)-(d’)] and cannot predict the electro¬ 
magnetic response for several layers [Fig. [5][d)]. Instead 
of isolated particles separated by the free-space gaps, 
we rather get a perforated high-index dielectric grating. 
Even though the period of the grating is smaller than 
the free-space wavelength, and thus no diffraction to the 
free space can happen, the guided mode in such dielectric 
membranes can be excited. The power coupled at some 
point is no longer localized within the particle, but can 
propagate within the dielectric layer, thus the metama¬ 
terial cannot be effectively homogenized. 

We have investigated the metamaterial made of high- 
dielectric rings. Even though the selected 5 = 100 
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is modest (dielectric constant of barium strontium ti- 
tanate may reach several hundreds [36|), it is represen¬ 
tative enough to make general enough conclusions. We 
may accurately anticipate the behaviour of the lower- 
index dielectrics. They might better suite to homoge¬ 
nization, having smaller refractive index contrast with 
the background. However, the unusual properties, based 
on strong electric and magnetic resonances, would also 
be weaker. 


In conclusion, the all-dielectric ring metamaterials can 
be described with electric and magnetic point dipole 
approximation if the meta-atoms size (ring radius) is 
smaller than 0.8 of the period. For the meta-atoms that 
almost touch each other, the dipole approximation is not 
applicable. We believe, the presented methodology and 
analysis results will be useful for development of low-loss 
all-dielectric based metamaterials. 
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FIG. 1. (Color online), (a) Dielectric ring metamaterial consists of rings arranged in a square lattice. The plane wave is incident 
perpendicularly to the layers. Electric field distribution in the unit cell for magnetic (b) and electric (c) dipole resonances excited 
with the incident plane wave. 
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FIG. 2. (Color online), (a) Absolute value of the reflection R 
for a single array of rings. Dimensionless electric (b) and 
magnetic (c) polarizabilities. The inner radius is fixed to 
Ri/d = 0.1. The outer radius R^jd— 0.2 (red), 0.3 (green), 
0.4 (blue) and 0.5 (purple). For clarity of presentation a verti¬ 
cal incremental offset is introduced to the spectra. The offsets 
of the spectra are marked with horizontal dotted line. This 
convention is kept throughout the rest of the article. The 
grey-shaded areas mark the range where higher-order modes 
are excited. 
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FIG. 3. (Color online), (a) Reflection R for a single array 
of rings and their dimensionless electric (b) and magnetic (c) 
polarizabilities. The outer radius is hxed R^jd — 0.4. The 
inner radius varies Ri/d = 0 (red), 0.1 (green), 0.2 (blue) and 
0.3 (purple). 
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FIG. 4. (Color online). Normalized electric (a)-(d) and magnetic (a’)-(d’) polarizabilities for R^/d — 0.2, 0.3, 0.4 and 0.5, 
respectively. The values of the actual period d* jd are (not on all graphs, only the cases of d* > 2i?2 are shown) 0.6 (red), 
0.8 (green), 1 (blue), 1.2 (purple) and 1.4 (orange). The real parts of the polarizabilities have positive offset for presentation 
clarity. The curves for d* jd — 0.6 (red) in (b) and d*/d = 0.8 (green) in (c) are multipled by 1/2 in order to fit the graph. 
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FIG. 5. (Color online). Absolute value of reflection \R\ for 2 
(red), 3 (green) and 4 (blue) layers of rings. Analytical predic¬ 
tions (thick solid lines) are compared to full-wave numerical 
simulations (thin lines with circles). The inner radius is fixed 
Ri/d = 0.1. The outer radius is R 2 /d = 0.2 (a), 0.3 (b), 0.4 
(c) and 0.5 (d). Grey area corresponds to the frequency range 
where high-order resonances are excited. 














































